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ABSTRACT 

We propose a new technique to directly measure the shapes of dark matter halos of galaxies using weak 
gravitational lensing. Extending the standard galaxy-galaxy lensing method, we show that the shape 
parameters of the mass distribution of foreground galaxies can be measured from the two-dimensional 
shear field derived from background galaxies. This enables the comparison of the ellipticity of the mass 
distribution with that of the light in galaxies, as well as an estimate of the degree of alignment between 
the stellar component and dark matter. We choose the specific case of an elliptical, isothermal profile and 
estimate the feasibility and significance of the detection of this signal. The prospects for applying this 
technique are excellent with large on-going surveys like the Sloan Digital Sky Survey. The expected signal 
is smaller, but comparable in significance to that of the mass in standard galaxy-galaxy lensing analyses. 
Since shapes of halos depend on the degree of dissipation and the transfer of angular momentum during 
galaxy assembly, constraints obtained from the analysis will provide an important input to models of 
galaxy formation. 

Subject headings: gravitational lensing, galaxies: fundamental parameters, halos, methods: numerical 



1. INTRODUCTION 

Current observations require the existence of dark mat- 
ter halos for galaxies. However, fundamental parameters 
such as their total mass and spatial extent are not well 
constrained. The mass distribution in galaxies is primar- 
ily probed via dynamical tracers of the galactic potential 
on various scales: stars in the inner regions, HI gas in 
regions outside the optical radius, and orbital motions of 
bound satellites in the outer-most regions (e.g. Zaritsky 
& White 1994). Probes of halo structure at radii devoid 
of any luminous tracers are therefore needed, weak gravi- 
tational lensing offers precisely that. 

Galaxy-galaxy lensing, the preferential tangential align- 
ment of the images of background galaxies around bright 
foreground ones, is detected statistically by stacking galax- 
ies. The first observational attempt to look for galaxy- 
galaxy lensing was made by Tyson et al. (1984). Recent 
studies have been very successful, and a signal at the 99.5% 
confidence level was first reported by Brainerd, Blandford 
& Smail (1996) using deep ground-based CCD data. Sev- 
eral subsequent studies using Hubble Space Telescope im- 
ages and ground-based data (Griffiths et al. 1996; Dell' 
Antonio & Tyson 1996; Hudson et al. 1998; Casertano 
1999; Ebbels et al. 2000; Fischer et al. 1999 [Sloan Digital 
Sky Survey (SDSS hereafter) commissioning data]) report 
unambiguous detection of a galaxy-galaxy lensing signal. 

In the method proposed in this letter, additional in- 
formation that is available but not exploited in current 
galaxy-galaxy lensing studies is utilized, namely, the light 
distribution of galaxies selected to be foreground lenses. 
The general results derived by Schneider & Bartelmann 
(1997) are used to relate the shear field to the mass mul- 
tipole moments. We show that the shapes and orienta- 
tions of the foreground galaxies (probes of the light) can 
be compared statistically to that of the shear field (probe 



of the mass), thus providing a direct method to compare 
the ellipticity of the light to that of the mass as well any 
potential misalignments between them reliably. These pa- 
rameters offer an important clue to the galaxy formation 
process, since they provide a quantitative measure of the 
importance of dissipation in the assembly of galaxies. Any 
variation of the flattening of the total mass (predominantly 
the dark matter component) with radius is a probe of the 
efficiency of angular momentum transfer to the dark halo 
and might provide insights into the structure and compo- 
sition of galaxy halos. With regard to the relative orien- 
tation of the light and mass in galaxies, the two compo- 
nents are likely to be aligned on average and any misalign- 
ments might occur transiently, for instance, after a violent 
merger. 

The outline of the paper is as follows: in the current 
status of modeling in galaxy-galaxy lensing studies is re- 
viewed. In §|| the formalism used to extract the shape of 
the mass distribution is presented, with the application to 
an elliptical isothermal mass model in §[I| The feasibility 
of detection of the flattening of the mass is examined in 
§5], and a measure of the alignment of mass and light is 
discussed in We conclude in §[7] with a discussion of 
the importance of studying shapes of dark halos for un- 
derstanding key issues in the formation and structure of 
galaxies and the relation between the luminous and dark 
component in galaxies. 

2. CURRENT STATUS OF GALAXY-GALAXY LENSING 

The galaxy-galaxy lensing signal is measured ideally 
from a deep image by selecting a population of brighter 
(assumed to be foreground) galaxies as lenses and mea- 
suring the induced shape distortion in the fainter (back- 
ground) galaxies. The shear signal obtained from di- 
rect averaging in radial bins around the bright foreground 
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galaxies is then stacked to obtain the radial profile of the 
shear 7 as a function of distance from the lens. This pro- 
vides reasonable constraints on the circular velocity of a 
fiducial halo, found to be in the range 210 km s _1 - 250 
km s _1 for a typical L* galaxy, but is found to be fairly in- 
sensitive to the radial extent of the halos (consistent with 
halo sizes in excess of 100h _1 kpc). In analyses to date, 
the errors are dominated by shot noise and error arising 
from the unknown redshift distribution of galaxies. How- 
ever, with forthcoming surveys like the SDSS, which plan 
to image over a hundred million galaxies in many bands 
and provide reliable photometric redshifts, the prospects 
for galaxy-galaxy lensing studies are extremely good (Fis- 
cher et al. 1999). 

3. EXTRACTING SHAPE PARAMETERS FOR THE MASS 
DISTRIBUTION 

The mass distribution of a lensing galaxy is described 
by the convergence field «(x), defined to be the surface 
mass-density E(x) expressed in units of the critical sur- 
face mass density E cr i t . The critical surface mass density 
depends on the precise configuration, i.e. on the angular 
diameter distances from the lens to the source D\ Sl ob- 
server to source D os and observer to lens D Q \. It is given 
by, E crit = Dls p ol ■ Standard galaxy-galaxy lensing 
provides a measure of the mass M within an aperture, 
which is given by 



(1) 



M= I drx w(x)k(x), 



where w(x) is a weight function normalized so that 
J d 2 x w(x) = 1. It is chosen to be continuous, differ- 
entiable and is required to fall-off rapidly to zero outside 
the aperture scale (3. The shape parameters of the mass 
distribution are characterized by the quadrupole moments 
of the convergence k(x) within the aperture (Schneider & 
Bartelmann 1997), which are defined as 



(2) 



Qij = / d x XiXj w(x) k(x), 



This tensor can be decomposed into a trace-free piece Q 
and a trace T defined as: 

Q = Q11 - Q22 + 2iQ 12 , T = Q n + Q 22 . (3) 

The ellipticity of the mass e K is then simply, 



e K = — 



T (al + bl) 



(4) 



where a K and b K are, respectively, the major and minor 
axes of the mass distribution, and (p K is its position angle 
relative to the positive x-axis. 

Schneider & Bartelmann (1997) have shown that the 
multipolc moments of k(x) computed from the observed 
shear _r(x) field. In particular, the quadrupole moments 
(Eq. ||) can be expressed as 

Q = f d 2 xe 2 ^ [g t (x) 7t (x) +ig x (x) 7x (x)], (5) 



where the rotated shear components 7* and j x correspond, 
respectively, to a tangential and curl-type shear pattern 



about the center of mass of the lens (see Rhodes, Refregier 
& Groth 2000 for an illustration) . They are related to the 
unrotated components by, 



7t = - [cos(2</?)7i + sin(2</?)7 2 ] 
7x = - [- sin(2^)7i + cos(2if)j 2 ] 



(6) 



where tp is the polar angle from the x-axis. The associated 
aperture functions gt(x) and g x (x) are given by 

g t {x) = 2V 2 {x) - x 2 w(x), g x {x) = -2V 2 (x), (7) 

where V n (x) = x~ 2 f Q dx' x' n+1 w(x'). Similarly, the trace 
part T and the mass M can also be written as 



d 2 xg t (x)~f t (x), M = / d 2 xht(x)-yt{x), (8) 



T 



where ht(x) — 2Vo(x) — w(x). 

4. APPLICATION TO THE ELLIPTICAL ISOTHERMAL 
MODEL 

We consider an isothermal model with concentric ellip- 
tical equipotentials (Natarajan & Kneib 1996). The pro- 
jected potential for this model is ip = ar, where a is the 
Einstein radius and r is a generalized elliptical radius. If 

the x-axis is aligned with the major axis of the potential, 

22 

the generalized radius is given by r 2 — + ; where e 
is the ellipticity of the equipotentials. The Einstein radius 
is related to the velocity dispersion of the galaxy a v by 
a = 47r(^) 2 (-^), and is of the order of 1" for galaxies. 

For a weakly elliptical model (e <C 1), the potential has 
the form, 



V>(x) 



a x\l 



cos 2(ip - ipo)}; 



(9) 



where tpo is the position angle of the potential, and reduces 
to that of a singular isothermal sphere in the circular limit 
(e = 0). Current observational limits on ellipticities of ha- 
los have been compiled in a comprehensive recent review 
by Sackett (1999), however, there are no constraints on 
the dispersion in the shape parameters, therefore, for the 
purposes of this calculation we have used the central value 
of e = 0.3. Besides, higher order terms will be approxi- 
mately (0.3) 2 = 0.09 which are still only corrections at the 
10% level, consistent with the assumption of small e. 

Restricting our analysis to the weak regime, to first or- 
der in e, the associated convergence k = V 2- 0/2, where 
ip = ar, is given by, 

«(x) - ^[1 + y cos 2(p- W) )], (10) 

and the complex shear 7 = 71 + 172 = [d 2 — d 2 + 2idi d 2 ]tp/2 
is 

7 - -^[1 + ycos2(^-^o)]e 2 ^, (11) 
The rotated shear components are thus 

It - ^[1 + yOob2( ¥ >- ¥ j )], 7x =0, (12) 
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yielding a tangential shear modulated by an elliptical pat- 
tern. 

The ellipticity of the underlying mass distribution k(x) 
needs to be related to that of the projected (2-d) potential 
</>(x). Using a normalized gaussian as the weight function 

w(x) = e~ x / 2/3 /(2-7T/3 2 ), we evaluate the integral for the 
quadrupole (Eq. |J) and monopole (Eq. E[) moments of 



M 



^-l^ r =V> (13) 



The ellipticity of the mass (Eq. ||) is thus e K = ^. The 
ellipticity of the potential e^, , obtained similarly by taking 
moments of ip, is = e/4. Note that the ellipticity of the 
potential computed above is smaller than e (by a factor 
of 4), as it is weighted by the circular gaussian window 
function. Comparing, the weighted ellipticities, we find 
that e K > e^, as expected since equi-potentials are always 
rounder than the mass contours. 

5. MEASURING THE FLATTENING OF THE MASS 
DISTRIBUTION 

We now show how these results can be used to measure 
the flattening of the mass distribution. As in ordinary 
galaxy-galaxy lensing, the galaxy catalog is separated into 
a foreground and a background subsample, using magni- 
tude, colors or photometric redshifts. The ellipticity of 
the galaxies in both subsamples is then measured in the 
by taking second moments of the light distribution. The 
ellipticities of the foreground sample yields the ellipticity 
of the light ei associated with each lens. While e; is ignored 
in ordinary galaxy-galaxy lensing, we instead align the 
foreground galaxies along their major axes before stack- 
ing. We then measure the average ellipticity of the mass 
e K as described above, by replacing the integrals in equa- 
tions (||) and (||) by sums over the sheared background 
galaxies. This yields a measurement of the component of 
the average ellipticity of the mass, e K y parallel to the that 
of the light, i.e. 

e„|| - Re(e:?,}, (14) 

where the ellipticities are taken to be complex numbers 
with e = ei + i&2, * denotes complex conjugation, and 
% — e //l e /| is the unit ellipticity of the light. 

We now compute the uncertainty in measuring e K n, by 
taking the square of the mean of the discrete estimators 
for M, T and Qn — Re(Q), and converting back into in- 
tegrals (Schneider & Bartelmann 1997). In the absence of 
lensing, we find 



a 2 [M] 



tibUfA 



a 2 [T} = 



n b n f 



A 



d 2 x h 2 (x), 
d 2 x g 2 (x), 
d 2 x [g 2 (x) + g 2 x (x)] , (15) 



where a 2 = (e 2 ) = (e 2 } is the variance of the intrinsic 
ellipticity distribution of galaxies (~ 0.3 2 ), rib and n/ are 
respectively the number density of background and fore- 
ground galaxies, and A is the area covered by the survey. 



For the elliptical isothermal model with the gaussian 
weight function, we can evaluate these integrals and find, 



a 2 [M] = 

° 2 [Q\\] = 



Airribn f Ap 2 

3<r 2 P 2 
AnnbUfA 



2o2 



(16) 



By propagating these errors in the definition of the ellip- 
ticity of the mass (Eq. 0), we can compute the signal to 
noise ratio (S/N) e = £ K \\/cr[£ K \\\ for measuring e re ii, and 
find, to first order in e, 



\(\ I I'-! i-, -r, i 1 i i -~ / 1 I \ I I llli 1,1,..,- 



5 arcmin 2 , 
3 /0.3\ / A 



(17) 



0.035arcmin- 2 / \a e J \1000deg 



Here, we have chosen to scale e K in units of 0.3, which is 
reasonable given current observational limits on the flat- 
tening of dark matter halos (see Table 3 in the Sackett 

(1999) review and references therein; Buote & Canizares 
1998). Additionally, in these scalings, we have used the 
survey specifications (ellipticity dispersion, number den- 
sity of foreground lenses nf and the number density of 
background galaxies n^, and approximate observed Ein- 
stein radius) quoted for the SDSS commissioning run pro- 
vided by Fischer et al. (1999), with a modestly expanded 
area (1000 deg 2 ) from the current area of 225 deg 2 . In 
the context of estimating the scatter arising in the mass 
estimates from galaxy-galaxy lensing due to halo shapes, 
Fischer et al. (2000) mention in passing that with 10 times 
more data than the commissioning run, halo shapes can 
be measured; this figure is comparable to our estimate of 
the S/N. 

Note that since these numbers have been taken from 
the SDSS commissioning run (which suffered from poor 
seeing), are conservative and do take into account several 
observational errors. The dispersion of 0.3 in the elliptic- 
ity distribution includes a correction for the seeing (see 
for instance weak lensing observations of Rhodes et al. 

(2000) and Bacon et al. (2000) and references therein). 
The effect of seeing and noise are also reflected in the 
modest number density assumed for background sources. 
Note, however, that in contrast to the case of ordinary 
galaxy-galaxy lensing, foreground galaxies will need to be 
sufficiently resolved to measure their shapes prior to align- 
ment and stacking. This could induce some degradation in 
the signal, but since the foreground galaxies are typically 
brighter and larger, this effect is expected to be small. 

The shape parameters of the mass will therefore be eas- 
ily detectable with SDSS in the near future (McKay et al. 
2000). For the total SDSS area of 10 4 deg 2 , the significance 
rises to 15<t. This in fact implies that potentially, even the 
radial dependence of the flattening can be studied by con- 
sidering annuli-shaped weight functions (for instance, the 
difference of two gaussians). Moreover, the degree of flat- 
tening as a function of the morphological galaxy type can 
also be studied. 

It is interesting to compare the (S/N) eK expected for 
measuring e K ii estimated above with that of the usual 
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galaxy-galaxy lensing (S/N)m — M/a[M) which measures 
the mass enclosed within an aperture. For the model con- 
sidered here, we find the following relation, 



0.17 



0.3 



(18) 



M 



Therefore, shape parameters can be measured with a sig- 
nificance which is smaller but comparable to that of the 
enclosed mass. Note that for the current SDSS survey 
area of A = 225 deg 2 , we find (S/N)m — 13 which agrees 
roughly with the significance of the reported Fischer et al. 
(1999) results, when averaged over all radial bins. 

6. MEASURING THE ALIGNMENT OF LIGHT AND MASS 

Given the good prospects expected from the above re- 
sults, one can be more ambitious and try to characterize 
the alignment of mass and light in more detail. We can 
make use of the amplitude of the light ellipticity e; , which 
we have not used until now. This can be done by group- 
ing the lens galaxies into several e;-bins, and computing 
e re ii separately for each bin. A more direct approach would 
be to consider the correlation function of the ellipticities 
of the mass and light, defined as 



C Kl = Re(e* K ei), 



(19) 



where the average is over all lens galaxies, and e K is an es- 
timate of the mass ellipticity of each lens derived from its 
associated background galaxies. While e K for an individ- 
ual lens is rather noisy, a significant measurement of C K i 
can be obtained by averaging over a large number of lens 
galaxies. This correlation function could also be computed 
for several annuli, and, would therefore, yield a direct mea- 
sure of the radial dependence of the alignment of mass and 
light. 

7. DISCUSSION 



The shapes of dark matter halos (see Sackett (1999) 
for a more comprehensive review) have been probed via 
many techniques and the consensus from these studies is 
that the precise shapes offer important clues to both the 
galaxy formation process and perhaps, even to the nature 
of dark matter. Cosmological N-body simulations suggest 
that dark matter halos are triaxial and that dissipation 
determines their shape. High resolution simulations find 
that the effect of dissipation (Katz & Gunn 1991; Dubin- 
ski 1994) is to transform an initially triaxial halo from 
prolate-triaxial to oblate-triaxial, while preserving the de- 
gree of flattening, yielding on average rounder and more 
oblate dark halos than those in dissipationless simulations. 

Comparing the shape of the mass profile inferred from 
X-ray data for a sample of ellipticals with that of the light, 
Buote & Canizares (1994; 1998) find that the dark matter 
is at least as flattened as the light and is definitely more ex- 
tended. The origin of the X-ray isophotal twist in the case 
of NGC720, they argue reflects an intrinsic misalignment 
of the stars with the dark matter. Keeton, Kochnek & Sel- 
jak (1997) incorporate the shape of the light distribution 
as a constraint in modeling individual lenses (that produce 
multiple images of background quasars) and find that an 
additional component of shear is required to match the 
observations. They speculate that this component could 
arise from a misalignment between the luminous galaxy 
and its dark matter halo. Our proposed technique will 
provide reliable measurements of the shape and orienta- 
tion of light and mass in galaxies, thereby aiding in the 
understanding of the coupling of baryons with the dark 
matter. 
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